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Photon propagation in a cold axion background with and without magnetic field 
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A cold relic axion condensate resulting from vacuum misalignment in the early universe oscillates 
with a frequency m, where m is the axion mass. We determine the properties of photons propagating 
in a simplified version of such a background where the sinusoidal variation is replaced by a square 
wave profile. We prove that previous results that indicated that charged particles moving fast in such 
a background radiate, originally derived assuming that all momenta involved were much larger than 
m, hold for long wavelengths too. We also analyze in detail how the introduction of a magnetic field 
changes the properties of photon propagation in such a medium. We briefly comment on possible 
astrophysical implications of these results. 

PACS numbers: 14.80. Va, 96.50.S-, 95.35.+d. 



I. INTRODUCTION 

Axions, originally introduced to solve the strong CP 
problem [ij], are to this date a viable candidate to consti- 
tute the dark matter of the universe 0. Their contribu- 
tion to the mass density results from the energy stored 
in the collective oscillations around the minimum of the 
axion potential 



a(i) = ao cos(mi), 



(1) 



with a frequency that is given by the axion mass to. We 
know that this mass must be somewhere in the range @ 



1 eV > to > l(T b eV. 



(2) 



The coupling of axions to photons takes place through 
the universal terrrQ 



C =n —-F F^ 

•'-077 — 1/077 „ „ r^ v r , 
27T } a 



(3) 



The constant <JVn is model dependent, but it is typ- 
ically of order 1 in most axion models [8(. The axion, 
being a pseudo Goldstone boson, satisfies the relation 
f a m ~ constant ~ fwTfi v , thus constraining the basic pa- 
rameters of the theory. However, the results presented 
below apply also to other light pseudoscalar particles, 
sometimes termed axion-like particles (ALP). The cou- 
pling between ALP and photons could in principle be 
stronger, since it is not related to their mass. 

Integrating by parts, we can write the term coupling 
axions or ALP to photons like 



1 



77 



with 



flu = Ti(t)S^o , 77 (t) = rj sin mi. 



(4) 



(5) 



The Lagrangian for a photon in the cold axion back- 
ground is then 



where F^ v = ^e^" a ^ F a p is the dual electromagnetic ten- 
sor. The dimensionful quantity /„ is the axion decay 
constant - the equivalent of f w as axions are assumed to 
be the pseudo Goldstone bosons associated to the break- 
ing of the Peccei-Quinn symmetry Upq (l)[l[. On f a we 
have a range of bounds: f a > 10 4 GeV coming from di- 
rect experimental searches of axions coupling directly to 
matterQ; /„ > 10 7 GeV from (somewhat weaker) as- 
trophysical constraints [(|, largely mass independent; or 
f a > 10 7 for 0.02 eV < to < 0.4 eV coming from the 
phase II of the CAST experiment Q. For some reviews 
of the experimental/observational search for axions see 

i- 



^ v F^<t> A , where 
are used intcr- 



1 This term is often written as Ca~, , 

(f>A is the axion field 4]. Both G^ 71 and g a -y-y 
changeably in the axion literature as coupling constants having 
dimensions E —1 . The constant g a u used here is howeever di- 
mensionless and it should not be confused with the latter. 



c = -\f^f^+ 1 -^a„f^, 



(6) 



and the relevant quantity to determine the physical effect 
of this coupling is 



Vo — 2ga~fj — 



a apTO 



7T fa 



(7) 



Now we can proceed to quantizing the photon field in 
such a background. This has been previously done in (Tl| 
in the case where rj(t) is assumed to be a constant, rj(t) = 
rjQ. It was found that in this case the two physical photon 
polarizations get their dispersion relations modified in 
the following way 



k 2 ±vo\k\. 



(8) 



As a consequence processes that are forbidden on 
Lorentz-invariance grounds such as 7 — > e + e~ or e — > 
have a non-vanishing probability if certain kinematical 
constraints are fulfilled. The interested reader can see 



2 



9] for possible observable consequences. If measured, 
these effects would constitute prima facie evidence that 
not only axions or ALP exist but they do constitute the 
primary ingredient of the dark matter of the universe. 

It was argued in d G3] that taking r)(t) as a constant 
was a good approximation if the momenta of all particles 
involved in the process were larger than m, the period of 
oscillations. However, if the wavelength of some of the 
particles are comparable or lower than the period of os- 
cillation one must necessarily deal with a time-dependent 
external potential. Thus it seems to us quite important 
to establish the basic principles of photon propagation in 
a time dependent axion background. For this reason in 
this paper we solve the problem of photon propagation in 
an oscillatory, but spatially constant, axion background 
exactly. We shall also include an external magnetic field 
to see how the combined effect modifies the properties of 
photons moving in such an environment. We will discuss 
at the end of the paper some possible physical conse- 
quences. 

To keep the paper technically simple we have approxi- 
mated the sinusoidal time dependence of the background 
by a square wave with the same period. A sinusoidal 
wave involves Mathieu special functions complicating the 
calculation enormously. We base this approximation on 
the similarity of the present effect with the emergence 
of the band structure in periodic potentials (l2|. exchang- 
ing time and space, and momenta and energies. It is 
well known in solid state physics that even such a simple 
model fully captures the esentials of metallic conductors 
and semi-conductors. Therefore we firmly believe that 
the physics of the problem being discussed remains un- 
altered by our technical simplification. 



II. SOLVING FOR THE EIGENMODES AND 
EIGENVALUES 

We introduce a Fourier transform with respect to the 
spatial coordinates only and write the photon field as 

( d 3 k 

A tl (t,S) = J _^ e *M M (t s fe). (9) 

The equation for A v (t, k) is 

'g» v {dl + P) - ie^T] a kp\ A v {t, k) = 0. (10) 
We now define 

S\=e^r la k p e llXpa T 1 Pk°, (11) 
which can also be written as 

=[(n ■ k) 2 - rfk 2 ] + fc 2 r/V 



-n 2 Wk v - (n ■ k){^k u + rfk»), 



(12) 
(13) 



and 



T —e^i la kp , S = S" M = 2 V 2 k 2 . (14) 
V 2o 



The properties of these quantities are discussed in [TT| j 
Note that the time dependence (due to r](t)) in P£" can- 
cels. With the help of these projectors we can write ([10] 
as 



A u (t,k) = 0. (15) 



g^{8 2 + k 2 ) + \j - (Pf - P^ 



To solve the equations of motion we introduce the polar- 
ization vectors defined in (Tl| and writcQ 



A v {t,k) = A(*MM). 

These vectors satisfy 

P^e lJ {k,±) = e^{k,±), P£%(£, T ) = 
and do not depend on t, so 

~d 2 + k 2 ± v m\]f±(t) = o. 



(16) 



(17) 



(18) 



As mentioned we will approximate the sine function in 
rj(t) by a square wave function: 



-T) 2nT <t<(2n + 1)T 



-770 (2n+ 1)T < t < 2nT 
The relevant parameters are 



^0 — 2g ail — 



a aom 



71" fa 



T = —. 

m 



(19) 



(20) 



There is an equation for each polarization. However, they 
are related. To recover one from the other we can just 
make the replacement 770 — > — 770- Also, because r)(t) 
changes sign after a time T in the square wave approxi- 
mation one solution is a time-shifted copy of the other: 
f-(t) = f+(t + T). In what follows we will work in the 
case A = +. It is obvious that the conclusions also apply 
to the other physical polarization, A = — . 

Since r](t) is piecewise-defined, we will solve the equa- 
tion in two regions: 
- Region 1: < t < T, rj(t) = rj 



d 2 fi(t) 
dt 2 



(k 2 + m \k\)f 1 (t)=0, 



f x {t) = A'e lat + Ae~ lat , a 2 = k 2 + m \k\. 
Region 2: -T < t < 0, r)(t) = 



d 2 f2{t) 
dt 2 



(k 2 - vo \k\)f 2 (t)=0, 



(21) 



(22) 



(23) 



2 When rj^ only has a temporal component, these polarization vec- 
tors actually reduce to the usual ones. 
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f 2 (t) = B'e i ^ t +Be 



i/3t 



P 2 = k 2 ~ Vo \k\. (24) 



We impose that both functions coincide at t — and we 
do the same for their derivatives 



/i(0)=/ 2 (0), /{(0) = ^(0). 



(25) 



We now write f(t) — e lut g(t) and demand that g{t) 
have the same periodicity as rj(t) 

gi (t) = e luJt h(t) = A'e^ u+a ^ + Ae l(ul ~ a)t , 
g 2 (t) = e iut f 2 (t) = B'e l ^ + ^ 1 + Be^'^, 

9l(T) = g 2 (-T), g[(T) = g' 2 (-T). (26) 

For these conditions to be fulfilled, the coefficients have 
to solve the linear system 

A' + A = B' + B 
aA! ~aA = f3B' - /3B 

(w + a)e l{u]+ ^ T A' + (w - a)e l ^ a ^ T A 

= {w + P)e-^+^ T B' + (u - P)e-^-^ T . 

(27) 

The linear system can be expressed as 



71/ 



( A '\ 








[')- 









(28) 



with 



M 1 



( 1 a 

1 —a 

-1 -j8 

V-l /? 



J(Lj+a)T 
i(u-a)T 



(u - a)e i{ ' 



(u-ot)T 

. e -i(a,+/3)T _( w + jS ) e -i(a, + /3)T 



. (29) 

The problem being discussed here is formally similar 
to the solution of the Kronig-Pennev|12j one-dimensional 
periodic potential, except the periodicity is now in time 
rather than in space. 

In order to find a non-trivial solution one has to de- 
mand the condition of vanishing determinant of M, which 
is 



cos(2wT) = cos(aT) cos(/3T) - 



/3 2 



2a(3 



sin(aT) sin(/3T), 
(30) 

with a and j3 given by (|22|) and ([24]) respectively. In order 
to get analytical expressions we will work in the limit of 
long wavelengths \k\T <C 1, which is just the one that is 
potentially problematic as discussed in the introduction. 
Expanding both sides: 



1 



irp2 



p. 



fc 4 



1 

12' 



which means 



12 



V 2 k 2 )T 2 + ..., (31) 



(32) 



If the determinant vanishes the system to solve is 

1(1 + |) 
h(a,/3,T) 



T e -i2u}T e ij3T 




B. 



.4 
B 

~B 



1 - 
1 



a — /3 e 



a + /3 e laT - e-^T e - l /3T' 



iaT e i2u)T _ & ipT 



(33) 
(34) 



iaT e i2uiT 



-ipT 
: YT e i2uT _ e i/3T 



a' a + (3 e iaT e i2 " T - e~^ T 



a~/3 e v 



a + /3 e iaT e l2 " T - e~^ T 



In the limit 770 < |fc|, \k\T < 1, 

B ~ B ~ 4|jfe|' B ~ 2|fc| 



(35) 



-^r. (36) 



Finally, imposing the usual normalization, 
/ h{t)fZ,{t) = 2K5{\k\-\k% 

we get 



(37) 



B 



k 2 + T]o\k\ 
2\k\+ m 



A 


2 


A' 




B 


+ 


~B 





k 2 - T] \k\ 
2\k\ - Vo 



B_ 
~B 



-1/2 



Vo 



(38) 



This completes the determination of the eigenvectors. 



A. Exact determination of the eigenvalues 



We can also solve (|30f exactly, without having to as- 
sume the long-wavelength limit as above, but this can 
be done only numerically. The solution only depends on 
770 and 771 through the dimensionless combination 71qT. 
There are values of k for which there is no solution, as 
seen in figure [1] However, these gaps get narrower when 
the product tjqT decreases. In practice, the largest pos- 
sible physical value for this quantity is tjqT = 1CP 14 and 
then the gaps are practically nonexistent and certainly 
totally irrelevant for the purposes of the present paper. 

It is interesting to investigate whether complex solu- 
tions exist for u in the forbidden narrow bands. We note 
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that the R.H.S. of the equation (|30|) is necessarily real, 
thus iv must necessarily be purely real or purely imagi- 
nary. In the latter case the L.H.S. is replaced by a cosh 
having as argument the imaginary part of 2loT. For this 
to have a solution, the R.H.S. must be positive and larger 



than one. Inspection of this term reveals that it is larger 
that one in the forbidden zones but actually alternates 
sign. Therefore not even an imaginary solution exists for 
the first, third,... forbidden regions. 




FIG. 1. Plots of the solution for rjoT = 1 and rjoT — 10~ . In the r/oT — > limit the solutions correspond 
to the straight lines ui ~ \k\ (plus their periodic repetitions). Small gaps develop but they become only 
physically significant when rjoT — 0(1). The physical region corresponds to the white area, the gray areas 
are just periodic repetitions. 



r 



B. Calculation of the transition e — ^ &y 

In order to make the photon field hermitian, we add © 
and its conjugate. Introducing creation and annihilation 
operators for each one of the proper modes we get (both 
polarizations are included) 



final state of an electron of momentum q and a photon 
of momentum k = p — q. 

(f\S\i) = iee* (fc, \)u ql »u p {2vf5^ (k + q-p) 

x / dtg*{t,k^ Ku+Eq ~ E ' }t (40) 



^3^[a(M)j(t,MK(M)e- b 



(27T) 

+a^(k,X)g*(t,k, X)e*Ak,\y kx 



(39) 



If we take rj(t) constant, g(t, k, A) = 1 and we have 
(f\S\i) =ieel{k,\)u q ru p {2K)H^ {k + q-p). (41) 



where kx = wt — k ■ x. Now we want to compute (f\S\i) In the square wave approximation (|19p . the time integra- 
tor an initial state of one electron of momentum p and a tion yields 

I 



</|S|i) = ieu^fupe*(k, \)(2n) 3 S{k + q- p)ir {AS(a + E q - E p ) 

+B5{(3 + E q - E p ) + A'S{-a + E q - E p ) - B'S(-f3 + E q - E p )} 



ieUq-fupe^k, A) I H — -=r (2^) J 5(fc + q - p)w i[ 1 - ) 5(a + E q - E p ) 



4|fc| 



2\k\ 



+5{f3 + E q ~E p ) + ^- i [6(-a + E q - E p ) + 6(-p + E q - E p )] 
A\k\ 



(42) 



(43) 



I 

Equation (I42[) holds for any value of \k\. The « sym- leading order in the 7/0 expansion this expression agrees 
bol indicates the use of (l36l). It turns out that at the 
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exactly with the one obtained in [1(| assuming that r]{t) 
was constant except for the fact that for each value of 
the polarization only one of the two delta functions that 
are not suppressed by terms of the form rjo/\ k\ can be si- 
multaneously satisfied; namely the one that implies that 



a or ft equals y fc 2 — \rj a \\k\, contributing with a factor 
1/2 with respect to what is found for constant r\ to the 
amplitude. Thus in the transition reduced matrix el- 
ement iA4 one gets for each polarization exactly one- 
half of what is obtained if r)(t) is constant. But in the 
present case both polarizations contribute so finally we 
get (1/2) 2 + (1/2) 2 = 1/2 of the result obtained with 
constant rj(t). 

As a consequence the predictions concerning the radi- 
ation yield of a high energy charged particle propagating 
in the cold axion background [9j are confirmed. 



III. PROPAGATION IN A MAGNETIC FIELD 



result can be obtained by direct inversion of the photon 
equation of motion (fTOj) . 

® ® ® 

i 

i i i 

= WW + 'WWW + 'WWW + ••• 

FIG. 3. Propagator in the axion background. 

We now compute the propagator in the presence of a 
magnetic field, using the second term in (1441) . In order 
to do that we use the propagator just found, represented 
by a double-wavy line and include the interactions with 
the external magnetic field. The dashed line corresponds 
to the axion propagator. Summing all the diagrams we 




We will now compute the propagator of the photon 
field with two backgrounds: a cold axion background and 
a constant magnetic field. To do so, we take ([3]) and write 
the axion and photon fields as a background term plus a 
dynamical field. We get two relevant terms 



C 



0,77 



±J» a %A v d a Ap + ??Zp?.ad l >A v F>»', (44) 

2 7T/o 



where a is the axion field, A^ is the photon field and 
pnv corresponds to a magnetic field: F " 1 — B l , i™ = 0. 
The first term is just ((4|). Here we will take ?y(i) to be 
constant; therefore the results that follow are valid only if 
the distance travelled by the photon, I, verifies I < 2n/m. 

The vertices and Feynman rules corresponding to these 
terms are shown in figure [2j With the first vertex we can 



® 



'WWVWV 




FIG. 2. The two relevant vertices. The corresponding Feyn- 
man rules are shown. 



FIG. 4. Full propagator after resummation of the interactions 
with the external B field. 



get 



where 



-ig- 



f = D F aX k\ 

Jp ^pa 1 "'A; 



m? + ig 2 K ' 



g-^P^, K = F^kpD^F^. 



(46) 

(47) 
(48) 



In order to simplify the result we shall assume that k-B = 
0, which may correspond to an experimentally relevant 
situation. Then we get 



fa = ik gl 



K = ik gl 



k 2 Bj - ir] {B x k), 
(k 2 - Vo \k\)(k 2 + m \k\) 



k 2 B 3 + irjQ(B x k)j 
(k 2 ~ V o\k\)(k 2 + Vo \k\) 



(49) 



(50) 



compute the propagator in an axion background, see fig- 
ure [3] The successive interactions with the axion back- 
ground can be summed up and the result is the propaga- 
tor 



= -i 



k 2 



P: 



(XV 



P 



/c 2 -77o|fc| fc 2 + ?7 |fc| 



(45) 

The physical polarizations, projected out by P± , exhibit 
poles at w 2 = fc 2 ± ryol^l as expected. The projectors are 



defined in (0 and = 



^\k\ 2 



Of course the same 



K = iklB 2 



k 2 



{k 2 -r,o\k\)(k 2 +r, \k\) 
and finally, defining b = gB, 

bjbi 



Vpu = Dp V + iklg^gl 



(51) 



(fc 4 - r] 2 k 2 )(k 2 - to 2 ) - k 2 k 2 b 2 



irjok 2 




x fc) 


-k(bx k)j 


- V 2 b 2 k 2 X 3l 


(fc 4 - 7? 2 fc 2 ) 


"(fc 4 


-r] 2 k 2 )(k 2 - 


to 2 ) — k 2 k 2 b 2 



(52) 
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A. Particular case: no axion background 

As a relevant particular case we now set ?yo = in the 
previous expression, i.e. we consider only the influence 
of the magnetic background, and get 



2V = D i*» + ik o9lgl 



b,b 



A k\k 2 -m 2 ) - k 2 k 2 b 2 ' 



(53) 



where now stands for the usual photon propagator, 
obtained after setting 770 = in (|4"5l) . 

This propagator has poles when fc 2 = k 2 and also for 



b 2 ±\ m 4 + b 4 + 2m 2 b 2 + Ab 2 k 2 



(54) 

If we assume that |6| is a small parameter and expand in 
powers of it, these poles in the frequency plane lie at 

fc 2 ~fc 2 (l + -^r ] + m 2 + b 2 , 



fc 2 , ~ k 2 1 



P 



(55) 



Physically this pole structure corresponds to the perpen- 
dicular polarization vector e± propagating unchanged, 
while the parallel polarization ey and the would-be lon- 
gitudinal polarization change their propagation^. 

For completeness we give the full propagator without 
the assumption k ■ B = 



^uy — Day -\- fnhv 



k 2 — m 2 + K 



(56) 



where D^ u = —ig^ v /k 2 is the usual photon propagator 
and 

fn = h n = -^(g^ob- k + g M Vk a ) (57) 



plane waves of well defined frequency w = fco and moving 
in the x direction. The Fourier transform with respect 
to the spatial component will describe the space evolu- 
tion of a photon state emitted at x — with polarization 
given by the vector €q. We decompose 



1 



A 



B 



C 



k 2 [k 2 (k 2 - m 2 ) - k 2 b 2 } k 2 -k 2 k 2 -F k 2 - G 

(60) 

where F and G are the roots of the denominator 



irr 



G = k 2 - — - - 
2 2 

2 1 



i 4 + Ak 2 b 2 



l-—\k 2 -m z 



F 



m 



4k 2 b 2 



h 2 



in agreement with (|55j). and 
1 1 



.4 
B 
C 



k 2 -Fk 2 -G 



1 



1 



k 2 -FF-G 



k 2 -GF-G 



k 2 b 2 



k 2 b 2 



3k 2 b 2 



(61) 

(62) 
(63) 
(64) 



Even for the largest magnetic fields conceivable the prod- 
uct b = gB is rather small compared to the range of 
acceptable values of the axion mass and it appears jus- 
tified to neglect C. The space Fourier transform of the 
propagator is then 



V^(k ,x) 



_ 9tw c ikt,\x\ _ fc ofr/A ^ c ifc,,|x| 



2k 2 



F 



(65) 



Let us now contract the propagator with the initial and 
final polarization vectors 



K 



1 

k 2 



{b ■ k) Z - b A k 



(58) 



Let us now restore the condition k ■ B = that is 
helpful in simplifying the formulae. In order to write 
the propagator in a more compact form we introduce a 
four- vector b^ — (0,6 = gB) 



ikob[jb u 



k 2 k 2 [k 2 (k 2 - m 2 ) - k 2 b 2 ] 



(59) 



Note the rather involved structure of the dispersion re- 
lation implied by (|53|). We consider the propagation of 



e^ u (k ,xy K e -^e ik °M + 

2ko 2ko 



(e- b)(e ■ b) ikM 



2k 



where b — b/\b\. Its squared modulus is 
\e»V^(k Q , X y Q \ 2 = 



Ak 2 



E 2 + 4 {ExE 2 + E 2 ) sin 2 



4m 2 



where 



(67) 



E-> 



(e- b)(e -b). 



(68) 



3 The labels X and || refer to directions perpendicular and parallel 
to the electric field, respectively, in the plane orthogonal to the 
propagation. 



This quantity, once properly normalized, describes the 
quantum mechanical probability of measuring the polar- 
ization represented by the vector e at a distance |x| from 
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the origin, where it was created with a polarization rep- 
resented by e . Since we restrict ourselves to the case 
k ■ B = and assume that the polarization vectors are 
orthogonal to the direction of propagation, we can write 

k = x, b = y. e = cosay+sinaz, cq = cos/3y+sin/3z, 

(69) 

so that 

Ei = cos(a — /?), E2 = cos a cos /3. (70) 
The extrema of (I57|) . for a given initial angle /3, are at 

[1 + 2/(2:)] sin 2/3 



tan 2a(x) 



4f(x) + [1 + 4/ (a;)] cos 2/3' 



where 



/(x) = sin 2 



4m- 



;fco|x| 



(71) 



(72) 



corresponding to the values of the angle where the prob- 
ability of finding the polarization vector is maximum or 
minimum. The mean value of the angle is 



a(x) 



1 



[l + 2/(ic)] sin 2/3 



2 [1 + 4/(2;)] +4/(x) cos 2/3' 



(73) 





If the electric field is initially parallel to the magnetic 
field, it remains parallel, i.e. a{x) = 0. Otherwise a 
rotation in the plane of polarization appears. 

The parameter characterizing the evolution is 
/co|2;|fe 2 /2m 2 . Usually fl3j mixing is treated via the clas- 
sical evolution equation 



(74) 



Note that the contribution from the Euler-Heisenberg la- 
grangian induced by the virtual contribution of electrons 
14] has not been included. It is not difficult to verify 
that both methods lead to the same dispersion relations 
in the case where 770 = 0. 



IV. INFLUENCE OF THE AXION 
BACKGROUND 

Now we return to the case with the axion background. 
In the limit where there is no magnetic field we recover 
the pole structure already discussed in the first section, 
u!± = k 2 ± 1 770 1 1 k I, whose implications in an astrophysical 
context were discussed in [9h11|. 

When the magnetic field is present, in addition to these 
poles we have three additional poles, manifest in the third 
term in (I52[) . that in the limit 7/0 ~ ^ correspond to 
the ones described following equation (j5"3")l . Taking into 
account that 770 is a very small quantity we shall disregard 



terms proportional to rj 2 in what follows. Then these 
three poles exactly agree with the ones discussed in the 
previous section. 

It is certainly of interest to consider how the simulta- 
neous presence of the magnetic field and the cold axion 
background influences the kinematics. If we consider the 
processes e —¥ ej or p — > Try discussed in the introduction, 
in the presence of both the cold axion background and 
the magnetic field the kinematical restrictions change. 
We defer this analysis to a separate publication. 

In order to consider the rotation of the polarization 
plane we note once again that 770 is a rather small pa- 
rameter. We shall neglect in the propagator all terms 
quadratic in 770; then (|52[) becomes 



2V(fc) 
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Now we contract with the polarization vectors: 
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where Ei and E2 are given in (|55| . (I70|) and 
E 3 = (e-b)[e -(bxk)]-(e -b)[e-(bxk)] =sin(a-/3) (77) 

We implement for the piece proportional to 770 a de- 
composition similar to the one described in the previous 
section. 
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F and G have been derived before. The new (tilded) 
coefficients are 
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We will again consider the propagation of an electromag- 
netic plane wave of frequency uj = kg in the x direction, 
perpendicular to the magnetic field. We have 

e"2V(fco,fciK = 
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Then 
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For small values of fcokl^/w 2 its square reduces to 
\e"V^(k ,x)e^\ 2 ^ 

P ^ 2 
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(85) 



[1 + 2f(x)] sin 2(3 + 3r) Q \x\ cos 2/3 

4/(x) + [1 + 4/(x)]cos2/3- 3770 |x| sin 2/3' 

(86) 



where terms of order 7/g have been neglected. 
The extrema of (1851) are at 



tan 2a(x) 



k \x\ , 



(87) 



16to 4 

and the mean value of the angle 

1 [! + 2/(x)]sin2^ + 3?7o|a:|cos2/3 
" ~ ~2 [1 + 4/(x)] +4/(x) cos 2/3 



V. PHYSICAL IMPLICATIONS 

In this paper we have seen that high-energy charged 
particles moving in a spatially constant but time varying 
axion background with frequency m radiate at a rate that 
agrees with the one computed in [t| [l(| , where the ap- 
proximation ui, k ^> m was assumed. It was seen in f9l. fioj] 
that the effect increases as the wave number of the emit- 
ted photons decreases, and its possible detection (if at all) 
is likely to occur in the MHz range of radiowaves. In this 
region k <C m and therefore the heuristic arguments used 
in [j| [l(| could be questioned. The calculation presented 
here settles the issue. The effect under discussion is quite 
real and to the best of our knowledge would constitute 
the clearest (perhaps even the onljo) observational evi- 
dence that axions or ALP constitute the bulk of the dark 



4 A possible exception would be the confirmation of 15] that ax- 
ions or ALP form Bose-Einstein condensates and caustics appear 
as a consequence in the galactic structures. 



matter component of the universe. The so-called direct 
observation experiments, such as ADMX, CAST or anal- 
ogous ones may find evidence for the existence of a parti- 
cle with the properties of the axion or ALP but this would 
not prove (although it would certainly be a tremendous 
hint) that axions or ALP constitute the missing mass of 
the universe. 

It is of course unfortunate that the amount of radia- 
tion predicted by the effect discussed here is very small; 
it is possibly within the sensitivity of long-wave radio 
antennae being commissioned or already deployed but 
around six orders of magnitude below the average value 
of the Galaxy synchrotron radiation background for the 
best value of 770. The effect approximately scales as 
Vo (^o/fc) 1 ' 5 - We expect this parameter to be ?/o < 1CT 20 
eV given the current bound for f a and the matter density 
(assumed to be due to cold axions). 

In principle observations in regions of low magnetic 
field could increase the signal/background ratio by sev- 
eral orders of magnitude as the synchrotron radiation is 
proportional to B 2 , assuming that the flux of cosmic rays 
stays at the average value in galactic regions of low mag- 
netic field. It should be noted that the assumption for 
the electron flux (electrons radiate most in the present 
mechanism |9|) was taken very conservatively to be the 
value measured by satellites, likely to be a gross underes- 
timate of the value in inner parts of the Galaxy. On the 
contrary, the background quoted is the observed value. 
In view of these considerations we believe that is impor- 
tant to refine the estimates before concluding whether 
this axion-induced Bremsstrahlung could be measurable 
or not, or used to place relevant bounds for 770 and hence 
on f a . Note that photons radiated via this effect are cir- 
cularly polarized, while synchrotron light is polarized in 
the plane of motion; measuring polarization may there- 
fore help in the detection of the effect. 

In any case it should be said that the effects under 
discussion could be considerably enhanced for ALP mod- 
els (assuming that the corresponding ALP particle con- 
densates similarly to axions proper) because some con- 
straints are evaded by these models. This is certainly 
something to have in mind and worth of further investi- 
gation. 

Another remarkable consequence of the presence of a 
slowly varying axion or ALP background is the fact that 
some wave-lengths (actually very narrow bands, see fig- 
ure I and its caption) are forbidden in the universe (or 
at least where there are substantial concentrations of 
cold dark matter, if this is constituted by cold axions 
or ALP). This opens of course a door for another line of 
experiments that could potentially probe these forbidden 
wave-lengths. The viability of these experiments, which 
appear very difficult unless the axion or ALP mass is 
known beforehand, deserves further investigation too. 

We have also studied the effect on the polarization 
of photons propagating in this oscillatory pseudoscalar 
background. We assume that rj(t) = rjo provides a good 
guidance. The results presented here have to be con- 
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sidered as exploratory and a more detailed account will 
be presented elsewhere. The relevant quantity that gov- 
erns the change in the plane of polarization is the ratio 
wxb 2 /2m 2 . The value of \b\ ranges from 10 -15 eV for 
magnetic fields of 10 T (such as the ones employed in 
CAST) to 10~ 6 eV for magnetar-strength fields, assum- 
ing that f a ~ 10 7 GeV. Taking m ~ 0.1 eV as a reference 
value for the axion mass this corresponds to the following 
approximate range 

10- 2 8 < A- < 10- 10 . (89) 
mr 

We qualitatively reproduce previous results [l3[ in the 
case where only the magnetic field is considered. How- 
ever, since we have derived the complete quantum propa- 
gator when photons propagate through an oscillating cold 
axion coherent background we can examine the modifi- 
cations due to it. We find that, quite remarkably, the 
modification in the ellipticity is independent of the light 
wavelength and also of the axion mass. It is probably 
even more notable that it is also independent of the mag- 
netic field itself, even if one needs to introduce one to 



begin with. We should warn the reader that because the 
results on the polarization are derived for constant r\ they 
are strictly valid for very short distances (< 2ir/m) only. 
The astrophysical consequences of the results presented 
here are yet to be fully explored. Note that to provide 
more realistic results other medium effects (such as an 
effective photon mass or the Euler-Heisenberg effective 
lagrangian) should be considered too. 

Clearly the presence of the cold axion background 
modifies the properties of photon propagation in rather 
interesting ways. The modifications are tiny, but some 
of them may perhaps be experimentally or observation- 
ally explored. This could possibly shed some light on the 
nature of dark matter. 



ACKNOWLEDGEMENTS 

We acknowledge the financial support MICINN 
through projects FPA2010-20807, Consolider CPAN as 
well as the support of SRG2009502. We thank A. Andri- 
anov and F. Mescia for numerous discussions. We would 
also like to thank K. Zioutas for encouraging us to exam- 
ine in detail the effect of magnetic fields. 



[1] R.D. Peccei, H.R. Quinn, Phys. Rev. Lett. 38 (1977) 

1440; S. Weinberg, Phys. Rev. Lett. 40 (1978) 223; F. 

Wilzcek, Phys. Rev. Lett. 40 (1978) 279. 
[2] L.Abbott and P. Sikivie, Phys. Lett. 120B, 133 (1983); M. 

Kuster, G. Raffelt and B. Beltran (eds), Axions: Theory, 

Cosmology and Experimental Searches, Lecture Notes in 

Physics 741 (2008). 
[3] R. Battesti et al., Lect. Notes Phys. 741, 199 (2008); P. 

Sikivie, Phys. Rev. Lett. 51, 1415 (1983) [Erratum ibid. 

52, 695 (1984)]; D. Dicus, E. Kolb, V. Teplitz and R. 

Wagoner, Phys. Rev. D 18, 1829 (1978); G. Raffelt, Phys. 

Rev. D 33, 897 (1986); D. Lazarus et al. Phys. Rev. Lett. 

69, 2333 (1992); Y. Inoue et al. Phys. Lett. B 668, 93 

(2008); G. Raffelt, Lect. Notes Phys. 741, 51 (2008). 
[4] K. Nakamura et al. (Particle Data Group), J. Phys. G 

37, 075021 (2010) 
[5] F. Bergsma et al. [CHARM Collaboration], Phys. Lett. 

B 157, 458 (1985) 
[6] J. Jaeckel, E. Masso, J. Redondo, A. Ringwald and F. 

Takahashi, Phys. Rev. D 75, 013004 (2007) 
[7] E. Arik et al [CAST collaboration], JCAP02(2009)008. 
[8] M. Dine, W. Fischler and M. Srednicki, Phys. Lett. B, 

104, 199 (1981); AR. Zhitnitsky, Sov. J. Nucl. Phys. 31, 



260 (1980); J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); 

M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. 

Phys. B 166, 493 (1980). 
[9] D. Espriu, F. Mescia, A. Renau, JCAP 1108:002 (2011). 
[10] A. Andrianov, D. Espriu, F. Mescia and A. Renau, Phys. 

Lett. B 684 (2010) 101. 
[11] A. A. Andrianov, D. Espriu, P. Giacconi and R. Soldati, 

JHEP 0909:057,2009; 
[12] R. de L. Kronig and W.G. Penney, Proc. Roy. Soc. Lon- 
don Ser. A 130, 499 (1931). 
[13] G. Raffelt and L. Stodolsky, Phys. Rev. D 37, 1237 

(1988). 

[14] SL. Adler, Ann. Phys. 67 (1971) 599; Z. Bialynicka- 
Birula and I. Bialynicki-Birula, Phys. Rev. D 2 (1970) 
2341; J.S. Heyl and L. Hernq uist, J. Phys. A 3 (1997) 
6485; W. Dittrich and H. Gies, hep-ph/9806417 H. Gies, 
|hep-ph/0010287] VA. De Lorenci, R. Klippert, M. Nov- 
ello, J.M. Salim, Phys. Lett. B 482 (2000) 134. 

[15] P. Sikivie, Phys. Lett. B567, 1 (2003). 



